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Stark effect is calculated by the perturbation theory method separately for the ortho and para 
water molecules. At room temperature, a 30%-difference in the energy change is found for the 
two species put in electric field. This implies a sorting of the ortho and para water molecules in 
non- uniform electric fields. The ortho /para water separation is suggested to occur in the course of 
steam sorption on a solid surface and of large-scale atmospheric processes. 

PACS numbers: 31.15.Md, 33.20.Ea, 33. 55. Be 
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It is a textbook knowledge that an isolated water 
molecule, H2O, exists in one of the two forms (spin iso- 
mers), ortho or para, with parallel or antiparallel proton 
spins, respectively In the first case a total magnetic 
moment of the molecule is 1, while in the second case 
it is 0. Under ambient conditions the gaseous ortho and 
para molecules are in a statistical 3:1 equilibrium (further 
referred as the normal O/P-ratio). In the gas, the ortho- 
para conversion due to collisions and radiation processes 
is highly improbable [2|. The valuable factor of conver- 
sion is magnetic interactions. In the absence of magnetic 
impurities, the estimated life-time of the gaseous ortho 
and para molecules is months Q ■ Thus, the water vapor 
is a true mixture of two gaseous ensembles consisting 
of magnetic [ortho) and nonmagnetic (para) molecules. 
This means that the ortho and para ensembles can in 
principle be spatially separated. 

For molecular hydrogen, the task of spatial separation 
of the spin isomers has been solved decades ago At 
present, the O/P separation of H 2 is usually attained by 
catalytically assisted cooling of gaseous H 2 at cryogenic 
temperatures. Separate existence of ortho and para wa- 
ter, however, is still questionable. The challenge is to 
pinpoint those physical properties of the chemically iden- 
tical ortho and para water molecules, which can be used 
for their separation. Here, by the perturbation theory 
calculations we study response of the ortho and para wa- 
ter molecules to the dc electric field (Stark effect), and 
find that these two species have noticeably different elec- 
trostatic properties. Although the study of Stark effect 
in water has a long history the effect, to our knowl- 

edge, has never been considered separately for the ortho 
and para modifications. 

As is generally accepted, we consider a water molecule 
as a rigid asymmetrical top and choose the coordinate 
system so that the £, ry, £ axes are directed along the three 
main inertial axes of the molecule [l[ . Further we follow a 
standard procedure described in details in Ref. [$]. The 
rotational part of Hamiltonian is written as: 



where J^,J v ,Jq are the angular momentum operators, 
and 1 a, Ib >Ic are the main inertial moments of the 
molecule. The symmetry of the Hamiltonian (TTJ) belongs 
to the D 2 group, thus the rotational energy levels are 
classified according to the irreducible representations of 
this group: A, B\, B2 and B3. After rotation by w along 
the C direction (CV(£)), the wave functions of the A and 
Bi symmetries do not change their signs, while the wave 
functions of the B2 and B3 symmetries do. 

In the Born-Oppenheimer approximation and neglect- 
ing the rotation-vibration interaction, the wave function 
of water molecule is |*) = \ip s ) ■ \ip e i) ■ \ip V ib r ) ■ \tp ro t), 
where \ij) a ), \ipei), \ipvibr)i an d iVVot) are the nuclear spin, 
electronic, vibrational, and rotational wave functions. 
According to the Pauli principle, the total ^-function 
should change its sign after permutation of the hydro- 
gen atoms. This permutation is equivalent to the C n (Q 
rotation. The functions \ipei) and \tp V ib r ) of the ground 
state are symmetric with respect to C W (Q. Whether or 
not \tprot) and \tp s ) change their signs, depends on in 
which states, ortho or para, the water molecule is. Obvi- 
ously, for the case of the parallel proton spins (the ortho 
molecule) the exchange of the protons does not reverse 
\ips). In this case, only |VVot) is responsible for the sign 
inversion of the total ^/-function. Therefore, \ip ro t) of the 
ortho molecule has either B2 or B3 symmetries. In con- 
trast, \ip s ) of the para molecule does change the sign when 
the protons (with antiparallel spins) are interchanged. In 
turn, |VVoi) remains sign-constant. Hence, the para rota- 
tional states have wave functions of the A and B\ types. 

We calculate the rotational levels and correspond- 
ing wave functions \ip r ot) using the perturbation the- 
ory method. The electric field is assumed to be weak 
enough so that the level shifts caused by the field are 
much smaller than the nominal distances between the 
levels. The direction of the z-axis of the immovable lab- 
oratory coordinate system (x, y, z) is chosen along the 
electric field vector. The interaction energy of the water 
molecule with electric field is: 
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V = -£-d z , (2) 
where d z is the projection of the dipole moment of the 
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water molecule to the z-axis. In the hrst approximation, 
the shifts of the levels are determined by the diagonal 
matrix elements of the operator ©. However, they are 
equal to zero for the asymmetrical rotators. Thus, the 
level splitting of the water molecule is a second order 
effect relative to the electric field: 
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where D = (J 1 ' Mjtn'\d z \JMjri) is a matrix element of 
d z , J is a total moment of the molecule, Mj is its projec- 
tion to the z-axis, the index n numerates different states 
with given J and Mj. The rotational wave function of 
the asymmetrical rotator is \ JMjn) = ^2 k C^ 7n) \JMjk}, 
where \ JMjk) are eigenfunctions of the operators J 2 , 
J z , and Jj, with eigenvalues, correspondingly, J(J + 1), 
Mj = J, J - I, -J, and k = J,J — I, - J. Then D 
can be expressed in terms of the eigen dipole moment of 
the molecule, S°\ directed along the £ axis, in a follow- 
ing way: 

D = Y, c k'' n ' ) * c k n) ( J ' M J' k '\^\ JM J k ) = 

kk' 

d {0) ■Y J C u' n> C { k 3n) {i) J - J '{-l) k - Mj ' x 

kk' 

v(2j + d(2J' + i).(_ j ; 

(-Mj Mj l)-*k*S*J*J» W 
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Klebsh-Gordon coefficients which are non-zero only for 
J' = J and J' = J ± 1 . In view of expression j4} , the 
general formula is reduced to: 
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(5) 

where aj n is the total shift of the "center of gravity" of 
the split level and (3j n is the shift of the sub-levels relative 
to this "center of gravity". From Eqs. and ©, by 
performing summation over Mj and Mji , we obtain the 
expressions for aj n and (3j n . The Stark behavior of the 
rotational levels of the water molecule is then calculated 
with the following parameters 0,11: a = h 2 /I A = 55.7 
cm" 1 , b = h 2 /I B = 18.6 cm" 1 , c = h 2 /I c = 29.0 cm" 1 , 
and = 1.84 Debye. 

The results of the calculations for J = 0; 1 are shown 
in Fig. [U The calculations are limited by the field of 
10 6 V/cm to meet the requirements of the second order 
approximation method we used. It is seen from the Fig- 
ure, that the singlet para-level \0A) with J = 0, shifts 
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FIG. 1: Stark effect for the water molecule: positions of the 
lowest four rotational energy levels, Ej, in external electric 
field. Light and bold lines show non-degenerated and twice 
degenerated levels, correspondingly. 



down when the field is applied. This is, of course, in ac- 
cordance with the general quantum mechanical theorem 
Q. The levels with J — 1 split into 2 sub-levels. One of 
these two sub-levels has Mj = and is not degenerated 
(thin lines in Fig. [T]). Another sub- levels correspond to 
Mj = ±1 and are twice degenerated (bold lines). 

Fig. [T] reveals qualitatively different trends for the or- 
tho and para levels in the field. When the field increases, 
the majority of the para- levels (the singlet \0A) and the 
twice-degenerated sub-level |l-Bi)) goes down. Only one 
of the para levels goes up. Thus, the net energy shift 
of the para levels is appreciably negative. In the case of 
ortho molecules, the twice-degenerated ortho levels (bold 
dash-dotted lines) diverge symmetrically as field goes up, 
the non-degenerated ortho levels lowering slightly only. 
As a result, the net energy shift of the ortho levels is 
notably, approximately two times, smaller than the net 
shift of the para levels. 

Taking into account the Boltzmann distribution of wa- 
ter molecules over the energy levels, one can write the net 
energy shifts AEo and AEp of the gaseous ortho and 
para molecules in the electric field £: 



= 1-2 E a Jn{2J + l)exp(-AEj Mj n/T) 

OP 2 E(2 J + l)ex P (-AE JMjn /T) ' 

(6) 

where AEjMjn — EjM jn ~ Eq are the energies of the or- 
tho and para levels Ej about the energy Eo of the lowest 
II-B3) and \0A) levels, correspondingly; T is the temper- 
ature; the summation is carried out separately over the 
ortho and para levels. 
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FIG. 2: Ratio of the energy Stark shifts for the para and ortho 
water molecules as a function of temperature. The dashed line 
notices the unity for eye. 



The ratio R — AEp / AEo gives a comparative re- 
sponse of the ortho and para water molecules to an ex- 
ternal electric field. R can be considered as a ratio of the 
forces acting on the para and ortho water molecules in a 
non-uniform electric field. As is seen from Eq. [6) it is a 
function of temperature T and does not depend on the 
field strength. 

We have calculated R(T) for 16 lowest rotational levels 
{J < 3), laying in the range of 0-283 cm^ 1 and contain- 
ing a 70 % deal of the molecules at room temperature. 
The Stark shifts of the higher levels of the asymmetric 
top decrease quickly when J increases therefore the 
contribution of the higher energy levels in R(T) is com- 
paratively small (less of 5% at room temperature). The 
result of our calculations is presented in Fig. O 

The drop of the R(T) at T < 30 K is easy to under- 
stand from simple considerations. Indeed, at low temper- 
atures the water molecules occupy mainly the two lowest 
ortho II-B3) and para \0A) levels. In accordance with Fig. 
[TJ the ortho | I-B3} level is of a superior "Stark flexibil- 
ity" . In an electric field it goes sharply down predomi- 
nantly causing the greater energy lowering of the ortho 
molecules. 

At elevated temperatures the additional para and ortho 
levels get filled with molecules resulting in levelling off of 
the total energy Stark shifts of the ortho and the para 
molecules. Surprisingly, however, the R(T) ratio docs 
not go to the unity, but exceeds it, reaching the 1.34 
value at room temperature. This implies that the para 
molecules in non-uniform electric fields experience a one 
third stronger electric force than the ortho molecules. 



The para molecules are thus proved to be more active as 
the "high field seeking" molecules. 

One can suggest the occurrence of the O/P sorting 
effect for the ortho and para water molecules in non- 
uniform electric fields, as it is typically happening with 
polar molecules 0]. Let us note, it is just electric rather 
than magnetic field that leads to a sizable distinction 
between the molecules with different magnetic arrange- 
ment. 

Obviously, one of the fields where the distinction be- 
tween electrostatic properties of the ortho and para wa- 
ter molecules can naturally be revealed is a water va- 
por sorption on a solid surface. The impinging water 
molecules are unavoidably exposed to the non-uniform 
surface electric field and, therefore, are suggested to un- 
dergo the O/P separation effect. A rough estimate [H of 
the total rotational energy drop A of the water molecule 
on the surface is: A ~ —d ■ £ ~ — 1 kJ/mol, where d 
= 1.84 Debye is a dipole moment of water molecule and 
£ ~ 10 6 V/cm is taken as an estimate of both, a typical 
surface electric field and the highest strength allowed by 
our perturbation method. As is seen, the dipole attrac- 
tion energy is not negligible in comparison with the real 
physadsorption energies of 5 - 20 kJ/mol 8]. The 30% 
difference in A for the ortho and para molecules looks 
large enough to be detectable during the adsorption. We 
suggest that this effect is responsible for the O/P wa- 
ter separation observed in absorption experiments with 
water vapor 

Beyond the sorption, a sizable effect of the O/P wa- 
ter separation can be suggested to occur in the atmo- 
sphere. In fact, the water vapor and atmospheric elec- 
tricity (ranging from 100 V/m to 1000 V/m normally and 
to millions V/m in lightnings) are inherent attributes of 
the atmosphere. Allowance of O /P water separation may 
give a new insight to the atmospheric phenomena. 

Summarizing, our calculations point to a remarkable 
difference in the forces acting on the ortho and para water 
molecules subjected to non-uniform electric field. This 
distinction is of a fundamental origin. Therefore, we sug- 
gest it should play an important role in nature environ- 
ments. 
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